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Abstract 

We consider the one-dimensional Cauchy problem for the Navier-Stokes equations with degen¬ 
erate viscosity coefficient in highly compressible regime. It corresponds to the compressible 
Navier-Stokes system with large Mach number equal to for e going to 0. When the 

initial velocity is related to the gradient of the initial density, the densities solving the com¬ 
pressible Navier-Stokes equations -pe converge to the unique solution to the porous medium 
equation [laiii]. For viscosity coefficient p{pe) = pf with a > 1, we obtain a rate of conver¬ 
gence of pe in L°°(0, T; for 1 < a < | the solution p^ converges in L°°(0, T; L^(M)). 

For compactly supported initial data, we prove that most of the mass corresponding to so¬ 
lution pe is located in the support of the solution to the porous medium equation. The mass 
outside this support is small in terms of e. 

1 Introduction 

The compressible Navier-Stokes equations in the multidimensional case with the Mach num¬ 
ber equal to read: 

' dtPe + div(p£'Ue) = 0, 

< dtiPsUe) + divipeUs (g) Ue) - div(2/r(p£)D(tte)) - V(A(pe)divri£) -F eVP{ge) = 0, (1) 

, &( 0 , x) = Qoix), PeUeiO, x) = mo{x), 

where Pe = Ps{t,x) and Ue = Ue{t,x) denote the unknown density and the velocity vector 
field, respectively, P{pe) = Qe, 7 > Ij denotes the pressure, D(ri£) = denotes 

the symmetric part of the gradient of u^, p{Pe) and \{pe) denote the two Lame viscosity 
coefficients satisfying 

p{pe) > 0, 2p{p^) + N\{p^) > 0, 

where N is the space dimension. 

Our purpose in this paper is to study the asymptotic behaviour of the global weak solutions 
to the one-dimensional Cauchy problem for (pQ), when e goes to 0, which corresponds to the 
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highly compressible limit. It is the opposite to the low Mach number limit leading to the 
incompressible Navier-Stokes equations. The latter has been intensively studied, especially 
in the case of constant viscosity coefficients. The first works in the framework of global weak 
solutions are due to P.-L. Lions and N. Masmoudi [25] and B. Desjardins and E. Grenier [S] 
(see also [lO]). The case of global strong solution with small initial data was considered by S. 
Klainerman and A. Majda in [2T] for well-prepared initial data and more recently in critical 
spaces for scaling of the equations by R. Danchin for ill-prepared initial data in [8]. Up to 
our knowledge the highly compressible regime has not been so well studied. The only results 
we are aware of concern Euler equations, see for example 0, m- The main difficulty to 
pass to the limit e —)• 0 in ([1]) with constant viscosity coefficients is due to lack of uniform 
estimates on {£»e}£>o (only the norm is conserved). In this case, it is not clear whether the 
sequence of weak solutions to ([T]) converges to solution of the pressureless system. However, 
in the case of viscosities satisfying particular algebraic relation 

\{Qe) = 2Qn'{Qe) - 2 h{q^), ( 2 ) 

additional compactness information is available. The first author proved in da HD that 
there exists a limit of the sequence {p^,Ue), called a quasi-solution, satisfying the pressureless 
system when N > 2. Moreover, the density of this system corresponds to the solution of 
the fast diffusion, the heat or the porous medium equation, depending on the choice of the 
viscosity coefficients. Relation ([2]) introduced by D. Bresch and B. Desjardins in m provides 
a new entropy structure which ensures estimate on the gradient of the density. In particular, 
it gives enough compactness information to pass to the limit when £ goes to 0. 

Stability of global weak solutions to © with ([2]) for s fixed was proved by A. Mellet and 
A. Vasseur [2^ using new energy estimate improving the integrability of the velocity. A 
particular choice of viscosity coefficients p{g) = g, X{g) = 0 satisfying ([2]) leads to the so- 
called shallow water system, for which the proof of existence of global weak solutions has 
been recently delivered by A. Vasseur and C. Yu [30], the extension to the case dJj) can be 
found in m, see also mm- Some ideas on construction of solution using the cold pressure 
(a pressure that is singular at the vacuum) can be found in |3l|29l|38] and in |1] using the 
notion of the K-entropy solutions. 

In the following paper, we extend nann by proving rate of convergence of g^ to a solution 
of the porous medium equation in a suitable functional framework. We restrict to the one¬ 
dimensional Euclidean space M, for which system ([T|) may be rewritten with a single general 
viscosity coefficient ^(0) > 0 in the form 


dtQe T dxige'^e') — 0 

dt (geUe) + dccigeuj) - dx{p{ge)dxUe) + £dxP{ge) = 0. 


Following m (see also m and 0 ), we rewrite the above system in terms of the effective 
velocity 

Ve = Ue + dxifige), where (p'{ge) = 

Qe 

which leads to a appearance of a parabolic term in the continuity equation 


dt Qe dx 


hiQe 

Qe 


dxQe ) T dx{g£V^') — 0, 


dt (QeVe) + dx(g£Ue:Ve) edxPige) = 0 . 


( 4 ) 
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This new regularising effect was observed in m in the framework of strong solution corre¬ 
sponding to initial density far away from the vacuum. The same concept was recently used 
in mu to construct regular solutions approximating global weak solutions to the full low 
Mach number limit system and to the compressible system. 

Our starting observation is that when ^^(O, •) = 0, the sequence {peVe}e>o converges to 0 
in a suitable functional framework. It heuristically implies that the limit p = lime_j.o Pe solves 
the fast diffusion, the heat or the porous medium equation 

I dtp - 9. = 0, 

I P(0,a;) = po{x), 

where 

p{p) = /5“ for a > 0. (6) 

The rate of convergence of {p^ — p) can be obtained by employing a duality method from [33], 
used by J. L. Vazquez to prove uniqueness of the very weak solutions to the porous medium 
equation. 

In the present paper we restrict to the case a > 1 because we are particularly interested 
in dealing with compactly supported initial data. When 0 < a < ^ the entropy discovered by 
D. Bresch and B. Desjardins (see m) allows to bound the density from below. Indeed, in m 
A. Mellet and A. Vasseur proved the existence of global strong solution to ([3]) for initial 
density far away from the vacuum. The main ingredient of the proof was to use entropy 
from m to estimate dx{p^ 2 ) in L°°(0,T;L^(M)) for all T > 0. This implies boundedness 
of p~^ in L°°(0, T; which allows to show sufficient regularity of solution in order to 

prove the uniqueness. This result has been recently extended by the first author to the case 
of the shallow water system [a = 1) [T5|. The main idea was to prove boundedness of 
appearing in (jH) by using the structure of the transport equation for It allows to bound 
p~^ by application of the maximum principle to the continuity equation. 

Finally, let us mention several important results concerning existence of global weak and 
strong solutions to system ([3|) with initial density close to vacuum. Presenting the exhaustive 
overview of the theory in this field is however beyond of the scope of this article. 

The existence of global weak solutions to system ([3]) was proven by Q. Jiu and Z. Xin 
in 120] . The main difficulty was to construct approximate regular solutions to (I3|) which 
verify all the entropies used in |33|. To construct such an approximation, Q. Jiu and Z. Xin 
employed the result on existence of global strong solution on bounded domain [—M, M] with 
Dirichlet boundary conditions. The latter result was proven by H-L. Li., J. Li and Z. Xin [23| 
who used energy estimate to control the Lipschitz norm of the velocity. 

There are also several interesting papers on system ([3|) with free boundary corresponding 
to the interface with the vacuum. The first result in this direction is due to D. Hoff and D. 
Serre m who proved the existence of global strong solution with discontinuous density at 
the interface for constant viscosity coefficients. An interesting extension of this result to the 
Neumann boundary conditions is due to P. B. Mucha |28j and for the reactive system see |23j . 
S. Jiang, Z. Xin and P. Zhang in m obtained similar results in the case of degenerate viscosity 
coefficients. The initial density considered in this paper is discontinuous at the interface and 
therefore is not compatible with the entropy introduced in |2| . We shall discuss this problem 
with more details later, in Section O 

The free boundary problem with initial density continuously connecting to the vacuum 
was analysed by T. Yang and H. Zhao in [36|. They proved the existence of weak solution 
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in a finite time interval in Lagrangian coordinates by showing that for sufficiently small 
T the L°°(0, T; L°°(n)) norm of p'^^^dxU is bounded. Unfortunately, this is not sufficient 
to come back to the Eulerian description since it is not clear how to solve the ordinary 
differential equation describing the evolution of the free boundary. In other words, the change 
of coordinates may not be well defined, see also [MIESIIST] . 

The article is structured in the following way. In Section [2] we formulate the notion of 
the weak solution to problem Q for e fixed and we state our main results. Next, in Section 
[3] we recall the result from [20] on existence of global in time weak solutions and derive the 
estimates independent of e. In the second part of this section we recall basic properties of 
the solution to the one-dimensional porous medium equation taken over from the monograph 
of J.-L. Vazquez |33] . In Section H] we prove the main results of the paper and in the end, 
in Section [5] we explain difficulties arising in the study of the free boundary problems and 
discuss possible extensions of our results. 

Notation: The letter C denotes generic constant, whose value may change from line to 
line. 


2 Main results 


Below we give the definition of a global weak solution to problem Q which we denote by 
(pejUg) or by (p, u) when no confusion can arise. The existence of such solution was proven 
by Q. Jiu and Z. Xin in [20|. We then formulate our main theorems in which we compare p^ 
with p the solution to the porous medium equation. 

Let us start with introducing the hypothesis on the initial data, we assume that 


Po > 0, mo = 0 a.e. on {x G M; po{x) = 0}, 

po GLi(M)nL°°(M), dxPo~'^eL^' 

|moP+'' 


(7) 


^ G L 

Po 




G 


where a > ^ and ly > 0 arbitrary small. 

Definition 2.1 Let e he fixed. We say that a pair of functions {p,u) is a weak solution to 
system ^ provided that: 

• The density p > 0 a.e., and the following regularity properties hold 

p G L°°(0, T; L\R) n L^{R)) n C([0, +oo), (iy^’°°(M))'), 

5^p“-^ G L°°{0,T;L\R)),^u G L^{0,T;L^R)), 
where (IU^’°°(M))' is the dual space o/IU^’°°(M). 

• For any t2>ti>0 and any tp G C^{[ti,t 2 ] x M), the continuity equation is satisfied in 
the following sense: 

f pfi{t 2 )dx- f pfi>{ti)dx= f f {pdtip + pudx'ip)dxdt. (8) 

i/]R. i/M i/lR. 

Moreover, for pv = pu+pdx<p{p), where ^'{p) = the following equality is satisfied 
f p'ip{t 2 )dx- f pip{ti)dx= f f {pdtfi + pvdxfi - pdx(p{p)dx'fi) dxdt. (9) 

i/M «.'M t\ t/M 


4 



• For any G C^{[0,T) x M) the momentum equation is satisfied in the following sense: 


/ mo^(O) dx + 

Jr 



{■\/p{y/pu)dtfi + {{y/puf + ep^)dxfi) dxdt 


-{p°'dxU,dxfi) = 0 , 


( 10 ) 


{p^'dxu, dxfi) 




y/pudxxjj dx dt 


2a 

2a - 1 




)y/pudx'fi dx dt. 


( 11 ) 


Below we give our first main result on the convergence of (pe, Ug) to a solution to the associated 
porous medium equation. 


Theorem 2.2 Let 7 > 1, a > 1. Moreover, assume that the initial data {pQ,mo) satisfy ^ 
and that 

= -V^dxifipo)- (12) 

1. System (0) admits a global weak solution {p^, ufij in the sense of Definition \2. h In addition, 
Pe converges strongly to p - the strong solution to the porous medium equation 


dtp - -dxxfi^ = 0 , 

a 

p{0,x) = po{x), 


(13) 


in the following sense, there exists a constant C > 0 depending on po such that 

\\{P-Ps){t)\\H-HR) <Ceh^. (14) 

2. For 1 < a < I there exists a constant C > 0 depending on po such that 

\\{pe - p)(i)||L 2 (iR) < Ceh^, (15) 


for all t >0. 


Remark 1 As mentioned above, the existence of weak solutions to system (0) was proven by 
Q. Jiu and Z. Xin in fMi- The novelty here is that our definition of weak solution incorporates 
additional equation ([9]). 

Remark 2 The initial data from m can be replaced by a condition which imposes this 
equality only approximately: 

< V (16) 

with T] > 0. It implies that the inequality 0 is true up to an error term coming from 

||(p-/>e)(t)||j^-i(R) < (17) 

for a constant C > 0 depending only on the initial data pQ. The same applies to inequality 

(03). 


mo 


\fpQ 


+ y/PQdxTiPo) 
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Remark 3 This result is restricted to the case 1 < a < | for technical reasons. This 
assumption allows to hound dxP and dxPe in L°°(0, T; L^(M)); actually we get even more that 
dxP^~^^‘^, dxP°‘~^^‘^ are bounded in L°°(0, T; L^(M)) for any a>\ (see also fT^). 

The next main result is the following. 

Theorem 2.3 Let 7>1, 1<Q!<|. Assume that the initial conditions satisfy ([7]) and that 
there exist two constants —oo < a < b < +oo such that 

supp[po] C [a,b]. 

Then, there exist constants C > 0 and —oo < oi < 6i < oo such that 

Qt ■= supp[p(t, •)] C [ai - Ct'^,bi + CT^]. (18) 

Moreover, there exists a constant C > 0 depending only on po such that 

\\p,it)ln^J\Li(^)<Ceh^l + t^), (19) 

where denotes the complement of Lit. 


Remark 4 The finite speed of propagation of the interface CHD is a classical result following 
from the theory of Porous Medium equation. 

Remark 5 Theorem \2.S\ provides a rate of convergence of {p^ — p) and so it extends the 
results of the first author proven in U3l\14^ when N > 2. Theorem 12.31 implies that the 
mass distributed outside Llf-the support of solution to the porous medium equation, is small 
of order £ 4 . The evolution of the support of Lit is known since the interfaces are described 
by the Darcy law (see Theorem VJ. 5\) . However, p^ may not be a solution to the free boundary 
problem. In other words, there is no reason for supppe^t, ■) to remain compact for all times 
t > 0. However, the mass that may be spread outside of Lit remains small. In this sense 
p{t,-) is a good approximation of p^{t,-) (for any t < max(l, i+^^) for c small) and 

we can think about the ’’quasi finite propagation” of the mass. 

3 Overview of known results 

This section is devoted to the overview and summary of results that will be used to prove 
our main theorems. 

3.1 Weak solutions to system (E]) 

We first recall the result on global in time existence of weak solutions to system ([3|) with e 
being hxed, whose proof can be found in the paper of Q. Jiu, Z. Xin ( [20], Theorem 2.1.) 

Theorem 3.1 (Existence of weak solutions) Let 7 > 1, a > 1/2 and let e > 0 be fixed. 
Assume that the initial conditions satisfy dZD- Then, system possesses a global in time 
weak solution {p, u) in the sense of Definition \2.1[ Moreover, this solution satisfies the fol¬ 
lowing inequalities uniformly with respect to £ 

p e C((0,r) X M), (20) 
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sup / pdx + max p < c, 

iG[0,T]iR (t,x)elO,T]xR 


( 21 ) 


sup f + (9a;(p“ 2 )H- dx 

t£[0,T]JR\ 7“1 / ^ 

+ [ f + A^) dxdt < c, 

Jo Jr ^ ' 

where c depends only on the initial data and A G L^((0,T) x M) is a function satisfying: 


Ap dx dt 


fl 


p" 2 y/pudx^p dx dt — 


2a 


2a — 


-/7 

1 Jo Jr 


(23) 


'9x(/o“ 2)^ti(^dxdt. 


Proof. The only element that we have to check is verification of the artificial formulation 
of the continuity equation Q. The rest was already proven by Q. Jiu, Z. Xin in [20] (see 
Theorem 2.1). Their proof combines the existence of global solutions on bounded domain 
[—M,M], proven in |24] with the diagonal procedure allowing to let M —)• oo. Therefore, we 
omit the details of construction of this approximate solution and focus only on justification 
that the formula (|9|) is valid. 

Let {p 5 ,us}s> 0 j Ps > C{5) be a sequence smooth approximate solutions to ([3|) with p{ps) 
replaced by psips) = piPs) + ^Ps^ ^ G (07/2), constructed in [2D] on the bounden interval 
Q = [—M, M] for M large, with the initial conditions 


P5(0) = Ps,o, Ps{d)us{0) = ms^o 


and the boundary condition 

UsldQ = 0. 

Assume that psfl, msfi converge to po, mo in the following sense 

P5,o^Po, strongly in L^(n), 

9x(7o 7 dx{Po ^) strongly in ^^(0) 


m 


5,0 


strongly in L^(fl), 
Po 


PS,o 

Wsfl 


2 + 1 / 


I mo 


2 + 1 / 


Ps ,0 


Po^'" 


G 


strongly in L^(fl). 


(24) 


We show that Q is weakly sequentially stable when (5 —>• 0. To this purpose let us recall 
the a-priori estimates derived in [23] following the strategy developed by A. Mellet and A. 
Vasseur in [26] for the multidimensional case. 

First recall that the classical energy balance gives for every T > 0 the following inequality 
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In addition, the BD entropy gives rise to the estimates 


+ 


P5- 

£ 




a — 1 


[ [ dxpjdxp'^ ^ dxdt 
Jo Jo. 

I 


+ 


< 


where we denoted p'^{ps) = ^ + 6p^^ Therefore, we obtain the following bounds 

^5 


£5 
e-i 

{us + dxP5{pi)f 


OJQ. 


dxpJdxPs dxdt 


(26) 


+ 


7-1 


p1 (0) dx. 


_ Q _^ 

\\y/p5Us\\L°°{0,T-,L^{n)) + \\dx{Ps ^)\\L^{0,T-,L^{n))+Hdx{Ps ^ ) IIl'^(0,T;L2(Q)) 

7 + q: —1 7+6 — 1 ^ 

+ 4P5\\1^^0,T ^L'lin)) d~ ^\\dxiPs ^ )llL2(0,r;L2(f2))+e5||6x(p5 ^ )\\L^{0,T-,L‘2(n)) ^ ' 

+ 11(^5^ + )dxUs\\L^(o^T;L^{n)) < C, 


where the constant is independent of 6. 
Following [23] we obtain 


L°°{o,T-,L°°{n)) < C. (28) 

This estimate can be then used to derive the improved estimate for the velocity, exactly as 
in [2^. Testing the momentum equation by us\us\'^ for u > 0 sufficiently small, and testing 
the continuity equation by '^'^ 2 +u summing the obtained expressions we obtain 


d f pslusl'^"^" 


-/ 

dt Jfi 


n 2 + 
= —£ 


dx + {u + l) / fis{P 5 )\u 5 \''\dxUs\'^ dx 


(29) 


' / dxpjuslusl" dx. 

Jq 

Due to the boundary conditions the r.h.s. can be integrated by parts and estimated as follow 

■ / dxPsUslusl" dx = £{u + 1) / p}\us\''dxus dx 
Jn Jn 


— £ 


< r 


e(z2 + 1) 


f ps\us\''\dxusf dx + f pp dx. 

Jn 2 Jq 


Note that due to our assumptions 27 — a > 0. For r sufficiently small, the first term can 
be absorbed by the l.h.s. of (|29p . while to control the second term, we can use the natural 
bounds for (l271) and (l28l) by choosing u such that ^ < 27 — a. Summarising, (f2^ yields 

sup \\ps\u5\‘^^''h^n) <C. 
tG[o,r] 

Having these estimates, we can follow the compactness arguments from [26] and from [23] to 
prove the convergence of the approximate solution {ps,us) to some weak solution {p^,Ue) as 
(5 —)• 0. Namely 


PS ^ Pe in C'([0,T] X D), 


[p's ^ j ^ j weakly in L^(0, T; l2(D)), 

V^U5 y/^Ue, PsUs Pe^e strongly in L^+5 (Q, T ; (O)), 

p'sdxUS ^ weakly in L^(0,T;L^(D)), 


(30) 
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where A satisfies (1231) . 

At this point, sequential stability of equations ([5]) and (fTU|l is verified. To justify the limit 
passage in ([9|), one has to prove the convergence of the term p^vs. Observe that we have 

P&V5 = P5U5 + psdx^sips) = PSU5 + ^ Ps^x (^Ps 1 

Due to (j3np the first term converges to PsU^ strongly in L^(0, T; L^(D)), while the second 
term converges to 2 a-i P^^^ weakly in L^(0,T;L^(D)); the last one converges to 0. 

In consequence 



P5V5 ^ PeVe weakly in L^(0, T; L^(D)), 
and so one can pass to the limit in ([9]). 

This is the final argument to prove the sequential stability of weak solutions to dH) on the 
bounded domain D = [—M, M] and with the Dirichlet boundary condition for u^. In order 
to let M Qo, we combine the diagonal procedure with the convergence of the initial data 
as it was done in [2D]. □ 


3.2 The porous medium equation 


In order to understand the qualitative properties of the limit solution to ([3|) , we recall several 
important features of the porous medium equations. The majority of them is taken from the 
excellent books of J. L. Vazquez [33], [32]. The porous medium equation can be written as 
follows: 

f dtp - dxxP'^ = 0, 

\ P(0,x) = po{x), 

with a > 1. 

In the sequel we shall set Q = (0, + 00 ) x M. Let us recall the notion of global strong 
solution for the porous medium equation (l3T|) (see [33] Chapter 9 for more details). 


Definition 3.2 ITe say that a function p G C([0, +00), 
to problem ^M\) if: 


positive is a strong solution 


p" G Lj^^{0,+oo-,L^{R)) and pt,dxxP°‘ S Ll^^{{0,+oo) x 
pt = pdxxP°‘ in distribution sense. 
p{t) —)• po as t ^ in 


The existence of global strong solution to (1311) is guaranteed by the following theorem (see 
page 197). 


Theorem 3.3 Let a > 1 For every non-negative function pQ G L^(M) nL°°(M) there exists a 
unique global strong solution p >0 of /T?!]) . Moreover, dtp G Lf^^{Q) for 1 < p < and: 


dtp > -- 


{a — l)t 
l|5tP(t)||Ll(R) < 


in V {Q), 


{a — l)t 
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Let Pi and p 2 be two strong solutions of (EIP in (0, T) x M then for every 0 < t < t 

||(P 1 - P 2 )_^(i)||Ll(K) < ll(Pl - P 2 )_^(t)||l 1 (r). (32) 

If Pi and p 2 are two strong solution with initial data (/Oi)o and {p 2 )o, such that {pi)o{x) < 
{p 2 )o{x) in M, then pi{t, x) < P 2 {t, x) for all {t, x) G ( 0 , +oo) x M. 

Likewise above, there exists a theory of global weak solutions with initial data being bounded 
measures. It covers the case of very important case of self similar solutions, the so called 
Barenblatt solutions U{t,x,M), of the form 

U{t,x,M) =t^F{xt^), F{C) = {C - , (33) 

where k = 2 iS^a+i) C > 0 with C = cM~ with 7 = 2 {a-i) ^ depends only on a > 1. 

The Barenblatt solutions verify ()31l) for t > 0 in the sense of distributions and with initial 
data po = Mdo, where do is the Dirac mass. 

Analogue solutions exist also in the case of fast diffusion equations corresponding to (l3Tl) 
with a < 1 . 

Remark 6 The comparison principle from Theorem \S.d\ ensures a finite speed of propaga¬ 
tion of solutions to the porous medium equation with compactly supported initial data, see for 
example WMi- In other words, the solution to the porous medium equation with compactly 
supported initial data remains compactly supported all along the time. Indeed, it suffices to 
compare such solution with the Barenblatt solutions. 


Let us now recall the so called — L°° smoothing effect (see [32] page 202). 

Theorem 3.4 Let pQ S L^(M). For every t > we have: 

p{t,x) < C'||Po|lil(o,T;R^)^ 

with a = /3 = and C > 0 depends only on a and N. The exponents are 

sharp. 

When po belongs also to L°°(R) the maximum principle holds. 


Let us denote by v the pressure in the sense of porous medium equation, i.e. 


V = 


a 

P , Vo = 


a — 1 




a — 1 


We recall an important theorem (see [33| page 376) which gives a precise behaviour of the 
interfaces when po has compact support. 

Theorem 3.5 Let x = s{t) denote the right interface of the solution p to (13ip emanating 
from the compactly supported initial data po G L^(M). For all t > 0 there exist the one side 
limits 

Dfv{t,s{t))= lim dxv{t,x), D~^s{t) = lim — [s(t + h) — s(t)]. (34) 

x^s{t)- h^o+ h 

Moreover, the Darcy law holds in the form 

= -Dfv{t,s{t)). 


The same result is true for the left interface. 


10 








4 Proof of main results 


Below we present the proofs of our main results stated in Theorem 12.21 and Theorem 12.31 We 
start by proving the rate of convergence from (1151) . the convergence of the approximate solu¬ 
tions (pe, Ue) to the solutions of the corresponding porous medium equation is a consequence 
of this estimate, the details can be found in [Ill- 


Proof of Theorem 12.21 

First let us note that thanks to Theorem 13.11 the continuity equation of system Q after the 
change of variables + dx(p{pe) with dx(p{pe) = dxPe reads 

Pe 

dtPe - -dxxPe + dxiPeVe) = 0, (35) 

a 

and it is satisfied in the sense of distributions. 

From (fT^ and (|26|) it also follows that 


sup \\y^Ve{t)\\L2(M) < ---T 

ie[0,r] (7-1)2 



< Cei 


(36) 


Next, we consider p - the solution to the corresponding porous medium equation with the 
same initial data po 


dtp - -dxxP'^ = 0 , 

a 

p(0,-) = po, 


(37) 


whose main properties were recalled in Section 13.21 
We now set 

Re = Pe- P- 

It follows from ([3^ and (l371l that Re satisfies the equation 

f dtRe - -dxxip"^ - p") + dxiPeVe) = 0, 

\ a 

y Re{^, x) = 0, X € M, 


(38) 


at least in the sense of distributions. 

Our goal is to estimate a relevant norm of Re in terms of e. To this purpose, we use 
duality technic in the spirit of J. L. Vazquez (see [33] Section 6.2.1). Testing (i38l) by -(/; G 
C'^((0, r] X M), we obtain 


[ [ [Redti^ + -(p" - p°')dxxi^)<^xdt 

+ f f {peVe)dx1p dx dt - f {Re'lp){T) 

J 0 J Jm. 


dx = 0. 


Let us now define 


a{t, x) 


Pe-P 


if Pe^ P 
if Pe = P- 


(39) 


Pe-P 

0 
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(40) 


This definition implies in particular that pf — p°‘ = aR^. We can hence rewrite (1391) as 
[ f {Reidt'tp +dxdt 

Jo Jr a 

+ [ [ {peVe)dr,'ll; dx dt - f (i2e'0)(r) dx = 0. 

J 0 ^ M. 1/ M 

The next step consists of solving the dual problem in the interval [—M, M] C M 

dtij +-andxx'J’= 0, (t,x) G [0,r] X 

a 

' ip{t,—M)=ip{t,M) = 0, tG[0,r], 

XT,x) = 0(x), xG(-M,M), 

where 6 G M)). Above, is a smooth approximation of a such that 0 < p < an < 

K < +00 (it will be precisely defined later on). This assumption guarantees that system (1411) 
is parabolic in the reverse time t' = T — t and it admits unique smooth solution if; on [0, T], 
Since ^ is a solution to (I4T]) . we can use (I4CT|1 to obtain 

[ Re{T)edx=- [ [ Re{a-an)dxxi^dxdt+ [ [ {peVe)dx'llj dx dt, 

Jr “ Jo Jr Jo Jr 


(41) 


therefore 


ReiT)edx 


< 


[ [ {PeVe)da 

Jo Jr 


■0 dx dt 


+ 


/ \Re\\{a - an)\\dxx'J’\dxdt. 

Jr 


(42) 


We now need to obtain the a priori estimate for dxxj’- To this purpose, we multiply equation 
(BH) by Cdxx'J’, where C = C(^) is a smooth and positive function such that ^ < C ^ 1 and 
dtC > c > 0, we get 


[ [ dtipCdxxip dx dt + - f [ (anidxxi^f dxdt = 0. 

Jo Jr “ Jo Jr 


(43) 


Integrating the first term by parts and using the fact that 0(0 •) G C^{[—M,M]) with all 
derivatives belonging to C'o([—M, M]) (see the general theory of [22], page 341), we obtain 

[ [ dti’Cdxxi’dx dt = - [ [ Cdxi’dxdti’dx dt =-^ [ [ dtidx'ilj)‘^Cdxdt, 

Jo Jr Jo Jr ^ Jo Jr 

{dxJjfdtCdxdt-^ [ {{dxJjfC){T)<ix + l- [ {{dxijfC){0)dx. 

E Jr Jr 


2 00 

Hence, it follows that 


[ [ dt'tpCdxx'iJ dx dt > ^ [ [ {dx'tp)‘^dtCdxdt-^ f {{dx'iljfC){T)(ix. (44) 

Jo Jr Jo Jr ^ Jr 


Plugging it into (|43]) we obtain 


ir [ {dx^pf 

^ Jo Jr 


1 


1 


dtC,dxdt-p- I {{dxO)'^C){T)dx + - I / Can(9xx0)^dxdt < 0. (45) 

J Jr ci Jo Jr 
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In particular, the assumptions on ^ imply that we have 


f [ {dxi>fdxdt +- [ [ an{dxx'ipfdxdt<c\\dx9\\l2(^). (46) 

Jo Jr ct Jo Jr 

Coming back to (H2|) and recalling that supp'0 and supp0 are included included in [0,T] x 
[—M, M] and [—M, M], respectively, we obtain 


dx 


[ Re{T)e 

Jr 

< C\\dx9\\i2(^^'^ 


rT rM 


^Ti 


'0 J—M 


(47) 


Re\ dxdtj + ||P£'ye||L2(o^T;I,2(K)) 


At this stage following ideas from [33] we construct the approximation an that is sufficiently 
regular and verifies ij < an < K for rj > 0 small and for K > r] large enough. This requires 
two steps of approximation. The first consists of taking aK,ri = min(Ar, max(7/, o)) with 
0 < r] < K, in other words 

{ K on {K < a{t, x)}, 
a on {rj < a{t,x) < K}, 
rj on {a{t, x) < rf\. 

Since aK,r) is T°°((0,T) x M), we know that ax^r^ is also in LP{{0,T) x (— M, M)) for 1 < 
p < + 00 . Therefore, in the next step, we can consider a smooth approximation a„ = 
o-K,r] * V'n (with ijjn a Standard regularising kernel) such that for n —)■ oo —)• aK,r} strongly 
in LP{{0,T) X (—M, M)) for all 1 < p < +oo. 

It remains to estimate the r.h.s. of (|171) . Since aK,rj > r] we deduce that also an > r] for 
any n, therefore 



2 rT rM 

^dxdt<— / \a — an\‘^\Re\'^ dx dt 

7 Jo J-m 


dr 


1 r dx dt + 


pT pM 

/ / ( 

JO J-M 


{{a-K)+ + r^f\R, 


2 


dx dt 


(48) 


We call the two integrals on the r.h.s. R and R, respectively. The integrand of R is pointwise 
bounded by 


a'^\Re\'^la>K + = (/?“ — Ps)'^da>K + v'^'i-a<K\Re\'^, 


thus, since (yo“ — p")^ is in L^( 0 , T; L^(M)) (which follows from (f 2 T]) and d?])), using the 
dominated convergence theorem we justify that X^(p“ — Pe )^la>K dxdt —)• 0 for K —>■ 
+ 00 . It implies that for K large enough, we have 



pT pM 

K)++ T])'^\Rirf dxdt < T]"^ + 1]“^ / 

Jo J-M 


\Rr 


< p2(i + 2TMC'), 


(49) 


where we used the fact that R^ is uniformly bounded in L°°((0,T) x M) for 0 < T < +oo. 
To estimate R, it suffices to observe that R^ is actually uniformly bounded in e in ((0, T) x 
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M) for any 1 < p' < +oo due to ([211) and due to L°° bound on the solution to the porous 
medium equation. We now use the fact that an —)• aK,r] when n —)• +oo strongly in ^^((0, T) x 
{—M,M)) to obtain that Jq |a — CLn\'^\Re\'^dx dt -^n^+oo 0. Using this, (P9]l and (H8l) 
we deduce that for n = n{r], K, M, T) large enough we have: 



d dji 


2 




Rs\‘^ dx dt < 2r]{2 + MTC). 


(50) 


Coming back to (I47p and using the Holder inequality we obtain that there exists C > 0 such 
that for every p > 0: 


From 


[ ReiT)9dx 

Jm 

< C'||(9a;0||L2(R) [ri{l + MT) + T 2 ||y^||ioo(o,T;L° 
and (|36l) we thus deduce that 


) 11 y/lhVe 11 (0,T;L2 (R))) • 


(51) 



< C'||93;^||L2(R)-||^==e2r2 ||■y/^||Loo(o,^;L°°(R)) 

< C\\dx0\\L2(^^'^£‘^T2 . 


(52) 


To conclude, let us observe that the above inequality holds for any 9 from C^{{—M,M)) 
with arbitrary M > 0 and that the constant C is independent of M. Thus, letting M to +oo 
we obtain 

WReiT^H-'^iR) ^ , (53) 

for any T > 0. 


Concerning the second part of Theorem 12.21 for 1 < a < |, since dxP^ ^ is uniformly 
bounded in L°°(0, T; L^(M)) and pi. is uniformly bounded in L°°(0, T; L°°(n)) we deduce 
that dxPe is uniformly bounded in L°°(0, T; L^(M)). Next, by ()53H and uniqueness of limits 
also dxP is bounded in L°° (0, T; (M)) (the same result can be obtained in dimension N >2] 
see [H]). Thus, by interpolation, the uniform boundedness of pe, p in L°°(0, T; L^(M)) and 
(l53l) . imply that 

\\R,{T)\\L2(^^<Ce^T-^. (54) 

This finishes the proof of Theorem 12.21 □ Note that for 1 < p < oo the rate of convergence 
in L°°(0, T; L^(]R)) can be proven by interpolation. 


Proof of Theorem 12.31 

Having proven the above rate of convergence, our next aim will be to provide some more 
information about the behaviour of p^ on the support of the solution to the corresponding 
porous medium equation. 

Under assumptions of Theorem 12.31 and from the properties of the solutions to the one¬ 
dimensional porous medium equation presented in Sectionwe deduce that the suppp(t, •) 
remains compact for all time t > 0 and it satisfies m- Indeed, since supppo is included in 
[a, b], we can consider the Barenblatt solution U{ti, M) at time > 0 such that 

I|7’o||l°°(R) < ||P(U,Tf)||j;^oo(R). 


14 











Last inequality is verified if M is large enough depending on ||/9o||l°° and ^ (see formula 
(1551) 1. Using the maximum principle we deduce that: 

1 1 

In particular suppp(t, •) is included in [oi — C'(t + ti)“+i, + U(t + ti)“+i] for some constants 

—oo < oi < 6i < oo and C > 0. 

Now, we may write the following sequence of equalities and inequalities 

IIP£(^)1o=IIl1(R) ^ IIpoIIlI(R) “ IIPeCOls^tllLhR) (fro™ tfro ™ass conservation) 

< l|p(i)l!^tllLi(R) “ IIPeC^ls^tllLLR) 

< ll(p(i) -/^£(i))lfitllLi(R) 

< |nt| 2 ||p(t) — Pe(t)||i 2 (]K) (from the Holder inequality) 

< Ceit^ (1 + t 2 (a+i) 'j ('from (fT8]l and ([551) 1. 

This finishes the proof of Theorem 12.31 □ 

5 Where are the interfaces? 

A natural question arising in the analysis of one-dimensional problems with compactly sup¬ 
ported initial data and free boundary conditions is propagation of the interface. From the 
classical theory for porous medium equation (see Section 13.21) we know that the interface 
moves with a finite speed. Moreover, according to ()19p we can estimate the mass correspond¬ 
ing to /?£ inside the support of p. This could suggest that the interfaces of free-boundary 
Navier-Stokes equations behave similarly to the interfaces of the porous medium equation. 
However, the issue of existence of global weak solutions to the free-boundary Navier-Stokes 
equations remains open. 

Below, we present an overview of partial results in this topic and explain the main prob¬ 
lems. 

1. A case of compactly supported initial density with jump discontinuity at the interface 
was studied in |19j . In this article S. Jiang, Z. Xin, and P. Zhang proved existence of global 
strong solution to ([5|) written in Lagrangian coordinates. Extending their solution outside the 
free-boundary domain by imposing the Rankine-Hugoniot condition on the velocity, allows 
to control dxU^ in L^(0, T; L°°(M)). It enables in particular to obtain global strong solution 
on the whole M. It is also sufficient to define the free boundary which verifies an ordinary dif¬ 
ferential equation. In this case one can pass from the Eulerian to the Lagrangian coordinates 
and conversely. 

It would be therefore interesting to compare the density pe from m with solution to the 
porous medium equation p emanating from the same initial data pQ. However, since pQ is 
discontinuous, it is not compatible with assumptions y/poOx^pipo) £ T^(M) and po £ L^(]R). 
As a consequence, we loose the a priori estimate on ^/p^Ve (1261) which was necessary to show 
damping of dx{peVe) in ([5]). In fact, the behaviour of pe is completely different than behaviour 
of Q. The discontinuity at the interface of p^ exists all along the time, while p becomes con¬ 
tinuous in arbitrary finite time (see Remark [6l). 
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2. The free-boundary problem with compact initial density, continuously connected to 
the vacuum was studied in [HB]. There, T. Yang, Z-a. Yao, and C. Zhu proved the local 
in time existence of weak solutions to Q in Lagrangian coordinates. Although this setting 
seems to be more adequate for our considerations, the existence of interface corresponding 
to this solution is not clear. Indeed, in order to determine the evolution of the interface one 
needs to come back to the Eulerian coordinates. Unfortunately, result obtained in [36] does 
not guarantee boundedness of the gradient of the velocity in L^(0, T; L°°(M)). This lack of 
control on dxU^ appears at the boundary, where the density vanishes. 

For the porous medium equation, the behaviour of the interface in the ID case is well 
understood since it verifies the Darcy law (see Theorem 13.Sp . It could be also appropriate 
to replace the usual free boundary condition for compressible Navier-Stokes equations by the 
Darcy law at the interfaces. Indeed in this case we could give a sense to Ue{t, si(t)) (where si 
denotes the right interface) by considering the limit lim 2 ,_,,^(^)- Ue{t,x) on the left hand side 
of the free boundary si. However with such condition the Lagrangian change of coordinates 
is a priori not possible. 
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